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We have studied the composition-induced metal-to-insulator transitions (MIT) of cation sub- 
stituted Lithium Titanate, in the forms Lii+^T^-^CU and LiAl y Ti2-j/04, utilising a quantum site 
percolation model, and we argue that such a model provides a very reliable representation of the non- 
interacting electrons in this material i/strong correlations are ignored. We then determine whether or 
not such a model of 3d 1 electrons moving on the Ti (corner-sharing tetrahedral) sublattice describes 
the observed MITs, with the critical concentration defined by the matching of the mobility edge 
and the chemical potential. Our analysis leads to quantitative predictions that are in disagreement 
with those measured experimentally. For example, experimentally for the LiAlj / Ti2- M 04 compound 
an Al concentration of y c ~ 0.33 produces a metal-to-insulator transition, whereas our analysis of 
a quantum site percolation model predicts y c » 0.83. One hypothesis that is consistent with these 
results is that since strong correlations are ignored in our quantum site percolation model, which 
includes the effects of configurational disorder only, such strong electronic correlations are both 
present and important. 
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The oxide spinel LiTi2 04 has been the subject of 
considerable experimental and theoretical study. It 
was first synthesised and structurally characterised in 
1971 by Deschanvres et alX Superconductivity, at 11K, 
was identified in 1973 by one of the present au- 
thors and his collaborators. 2 A comprehensive study 
of the normal state and superconducting properties of 
Lii +2; Ti2- x 04 (for < x < 1/3) was reported in 1976^4 
and a superconducting transition temperature of 13K was 
observed. A recent review 5 highlights many of the ad- 
vances made since then. 

There are several reasons to study this system. Firstly, 
it is interesting to note that superconductivity among 
spinel systems is very rare; e.g., of the 300 or so 
known spinels^ only four of them are superconductors 
- CuRh 2 Se 4 (T c = 3.49 K), CuV 2 S 4 (T c = 4.45 K), 
CuRh 2 S 4 (T c = 4.8 K), and LiTi 2 4 (T c = 11.3 K) - 
and only one of these four is an oxide. So, that oxide, 
LiTi 2 04, has the highest transition temperature of any 
spinel. Secondly, conduction in this system is believed 
to take place on the Ti sublattice via the i 2g orbitals, 
as suggested, e.g., by electronic structure calculations^ 
and these sites form a corner-sharing tetrahedral lat- 
tice(CSTL). Thus, this system represents an example of 
conduction on a, fully frustrated three-dimensional lattice. 
Also, in this paper we will argue, supported by consider- 
able experimental evidence, that the conduction paths of 
Lii +2; Ti 2 _ I 04 and LiAl a Ti 2 _ y 4 are excellent physical 
realizations of quantum site percolation^^ 

Furthermore, and central to the motivation for our 
work, these same electronic structure calculations 6,7 
point out that this is a narrow band electronic system, 
with the bandwidth of the i 2g bands of the order of 2-3 
eV, thus suggesting that perhaps strong electronic corre- 



lations are present and important. Indeed, others have 
reached similar conclusions; notably, the phase diagram 
of Alex MiillerjAfi summarising a view of how the in- 
creased strength of electronic correlations in transition 
metal oxides leads to higher and higher superconduct- 
ing temperatures, includes the Lithium Titanate sys- 
tem. Although the original experiments and analysis sug- 
gested a weak-coupling BCS-like s-wave superconductor 
it was later suggested 11 that off stoichiometry this mate- 
rial is in fact an "anomalous" superconductor (although 
this claim is not without criticism 12 ). We also mention 
that photoemission studies of Edwards, et alM are in- 
terpreted in terms of strong correlations, and magnetic 
susceptibility^ and specific heat dataii are interpreted in 
terms of a density of states that is moderately to strongly 
enhanced (see Refii 5 - for a discussion of these and other 
experiments). Taken together, these experimental results 
form a reasonably strong case for the presence of elec- 
tronic correlations that are important to the physics of 
these materials. 

Lastly, we mention the recent discovery of the first d- 
electron heavy fermion compound, LiV 2 04^ This sys- 
tem also assumes the spinel structure, but so far no su- 
perconductivity has been observed. The active transition 
metal ion in LiV 2 04 has a formal valence of d 1 ' 5 , whereas 
for LiTi 2 4 one considers d 5 ions. Thus, LiTi 2 04 is a 
lower electronic density system than LiV 2 04, and an un- 
derstanding of its behaviour would seem to be a prereq- 
uisite to a full understanding of the Vanadate material. 
For example, why does LiTi 2 4 superconduct, with a rel- 
atively high T c , while LiV 2 4 does not superconduct at 
all? 

In attempt to gain more understanding of the 
LiTi 2 04 system, and, in particular, to try and un- 
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derstand whether or not strong electronic correla- 
tions are present, we have examined the density 
driven metal-to-insulator transition (MIT) of the re- 
lated Li 1+2 .Ti2_2;04 and LiAlj / Ti 2 _ l/ 04 compounds; for 
xmit ~ 0.12 and Dm it ~ 1/3, transitionoVMVV 9 to 
a non-metallic state (which we refer to as insulating) are 
encountered. To be specific, we use a one electron ap- 
proach to study this transition employing a quantum site 
percolation model. Our work may be viewed as address- 
ing the question of whether or not the MIT undergone 
by this system is driven by disorder only, similar to an 
Anderson-like MIT. We find that the answer is no, and 
thus this work provides indirect theoretical support for 
the proposal that strong electronic correlations are im- 
portant in a description of the complicated transitions 
undergone in the LiTi2 04 class of materials. 

To fully explain our model we note the following: (i) 
Electronic structure calculations^ show that the bands 
arising from the Ti 3d orbitals are separated from the 
O 2p band by about 2.4 eV; thus, the electronic va- 
lence state may be represented as Li +1 (Ti +3 - 5 ) 2 (0 _2 )4, 
and we ignore the oxygen sites and focus on only the 
Ti sites. The crystal octahedral field around Ti cations 
splits the Ti 3d bands into two separate and nonoverlap- 
ping t2 g and e g bands, with the e g bands split off above 
the t2 g bands. Thus, formally this is a very low filling 
system — l/12th filling of each of the (approximately) 
degenerate t2 g bands. Although we have generalized our 
work to include all three ti g bands, here we will present 
results for a one-band model of the Ti sublattice, and 
thus the stoichiometric compound is represented by a 
1/4-filled band, (ii) Crystallographic refinements of the 
excess Li system Lii +:c Ti2_ a; 04 and the doped Al sys- 
tem LiAlj,Ti2_j,04 have consistently demonstrated that 
both the excess Li and doped Al ions enters substitution- 
ally onto the Ti sublattice (octahedral sites of the spinel 
structure) ii^*i2i22i2i Assuming that the Li/Al ions that 
are substituted into the corner-sharing tetrahedral lattice 
are fully ionised, these sites would block any conduction 
electrons from hopping onto such sites; e.g., a simple ar- 
gument supporting this follows from noting that the Li + 
ions will be at least doubly negatively charged relative to 
the occupied Ti 3+ and unoccupied Ti 4+ sites that would 
exist in the absence of substituting Li, and thus electrons 
will avoid these sites in favour of the Ti sites. We shall 
assume that these Li-substituting sites are removed from 
the sites available to the conduction electrons, which im- 
plies that this system represents an excellent physical 
realization of site percolation. 

Using such a model the simplest approach to character- 
ising the MIT would be to then associate the transition 
with the critical concentration at which the classical per- 
colation threshold is reached. For corner-sharing tetra- 
hedral lattices, we have completed a large scale Hoshen- 
Kopelman 2 ^ search, and have determined that this con- 
centration corresponds to a probability of finding an oc- 
cupied site at the transition of p c ~ 0.39 ± 0.01. Noting 
that the relationship between the probability p of site 



being occupied by a Ti ion (in the stoichiometric Ti sub- 
lattice), and the excess Li concentration x, is given by 
x = 2(1 — p), with an identical y = 2(1 — p) relation 
for Al added to the Ti sublattice, in contrast to previous 
statementsiii^ this p c corresponds to a critical excess 
Li, or added Al, of x c = y c = 2(1 — p c ) ~ 1.2; that is 
Li2.2Tio.8O4 or LiAl1.2Tio.8O4. These very high levels 
of doping are well beyond the observed Xmit ~ 0.15 or 
Umit ~ 0.33 concentrations at which the MITs occur. 
In fact, such a system would require large positive Ti va- 
lencies well beyond anything seen in nature! Thus, the 
physics of these MITs is more complicated than simply 
the loss of an infinite maximally connected cluster at p c . 

We now consider a more accurate model for this sys- 
tem, a so-called quantum site percolation model, which 
includes the dynamics of the electrons hopping on the 
conducting, disordered sublattice of Ti sites. To be spe- 
cific, the near-neighbour tight-binding Hamiltonian for 
LiTi2 04 is 

w = E e ^ - 1 E ( c ta + h - c -) (!) 

i <ij> 

where i labels the sites of the (ordered) Ti sublattice, cj 
(cj) is the creation (annihilation) operator of an electron 
at site i,J2<ij> represents the sum over all nearest neigh- 
bour sites of a corner-sharing tetrahedral lattice, and t 
is the near-neighbour hopping energy. To produce our 
model of quantum site percolation for the doped systems, 
the on-site energy Ej is determined by the probability of 
occupation, denoted by p, of a site being either a Ti ion, 
or a Li or Al dopant ion: 

P(s 4 ) = p5{ei - e T i) + (1 - P)5{e t - ex) (2) 

where exi (£x) is the on-site energy when an electron oc- 
cupies a Ti (X = Li or Al) site. In order to enforce that 
itinerant electrons move only on Ti sites we set £xi = 
and £x —> 00, and this limit connects this system with a 
quantum site percolation Hamiltonian. 8-9 Such consider- 
ations lead to the introduction of the quantum percola- 
tion threshold, 9 usually denoted by p q . To be specific, p q 
is reached when all single-electron energy eigenstates of 
the above tight-binding Hamiltonian are localized. That 
means that p q is always larger than p c since the presence 
of extended states necessarily requires an infinite maxi- 
mally connected cluster. Further, the evaluation of this 
quantity is warranted since p q > p c implies that the the- 
oretical predictions of x c and y c given above would be 
reduced by quantum percolation. 

Our evaluation of p q for the corner-sharing tetrahedral 
lattice proceeds as follows. We have considered various 
realizations of site percolated lattices for several system 
sizes for a range of dopant concentrations; to be spe- 
cific, we consider lattices of size two, four, six, and eight 
cubed conventional unit cells (noting that there are 16 
Ti sites in the ordered lattice per conventional unit cell), 
and then for each p we examine 100, 50, 20 and 10 re- 
alizations consistent with this p, for two, four, six, and 
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eight cubed lattices, respectively. For each realization we 
first apply the Hoshen-Kopelman algorithm to identify 
the maximally connected cluster. Then, we diagonalizc 
the one-electron Hamiltonian describing the electron dy- 
namics on this cluster. Then, to determine the localized 
vs. delocalized nature of the single-electron wave func- 
tions for the maximally connected cluster, we have used 
the scaling of the relative localization length as a func- 
tion of system size, as described by Sigeti et alm& This 
localization length, for a particular eigenstate, is defined 
by 



(3) 



where \ipi\ 2 is the probability amplitude for this eigen- 
state at site i, and d(i,j) is the Euclidean distance be- 
tween lattice sites i and j. Then, the relative localization 
length is just the ratio of this quantity to that for a state 
having a uniform probability amplitude throughout the 
entire maximally connected cluster (we denote the latter 
by Ao) — this ratio thus provides a useful measure of the 
effective size, or localization, of a particular eigenstate 
relative to a Bloch state on the same maximally con- 
nected cluster. The utility of this quantity (and we will 
describe its use for another problem below) is that if the 
quantity decreases as the system size is increased, that 
eigenstate corresponds to one that is spatially localized; 
the opposite behaviour is expected for delocalized states. 
We have used this quantity as a means of identifying p q . 

As a test of this method, we note that for three- 
dimensional lattices, reliable estimates exist only for the 
simple cubic lattice, and these were obtained with a va- 
riety of different methods — e.g., see the discussion in 
Ref^i A value of p q = 0.44 ± 0.02 was identified^ and 
we have found that our method reproduces this number. 

Using this method we find a value of p q for corner- 
sharing tetrahedral lattices (with near-neighbour hop- 
ping only) of p q = 0.52 ± 0.02, and if we then as- 
sociate this quantity with the concentration at which 
the MIT occurs in doped LiTi 2 C>4, one finds x c = y c = 
2(1 — p q ) = 0.96, which correspond to Li1.96Ti1.04O4 and 
LiAlo.g6Ti1.04O4. Again, these concentrations are much 
higher than the experimentally measured values. 

The reason for both of these failures is clear and has 
been suggested beforeii — as the concentration of doped 
cations is increased, the density of available 3c? electrons 
is decreased. In fact, for xbi = 0.33 and ysi = 1-0 these 
materials become insulators simply because the bands of 
all allowed states are empty (we refer to such a state 
as a Band Insulator (BI)). So, if a one-electron descrip- 
tion is going to correctly reproduce the MIT, we must 
account for the changing of the electronic density with 
cation dopant concentration. (This notwithstanding, we 
needed to determine what values of x c and y c were pre- 
dicted by p q in case they were less than either xbi, or 
Vbi, respectively. Clearly, they are not.) 

Also, since the above arguments point to Ti occupa- 
tion probabilities well above the quantum percolation 
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FIG. 1: A sequence of schematic diagrams that summarize 
how the MIT in Lii+ a; Ti2-2;04WOuld proceed if the transition 
was caused by disorder only. The densities of states vs. en- 
ergy are shown for several x. The hatched regions identify the 
location of extended states, and the unhatched (solid) regions 
of density of states represent the location of localized states, 
(a) In LiTi2C>4all eigenstates are extended and the system is 
a i-filled d band conductor, (b) By doping Li cations ran- 
domly into the Ti sites the system becomes disordered; thus, 
some energy eigenstates near the band edges become local- 
ized, and, in particular, all states below the mobility edge 
\im are localized. At the same time the density of itinerant 
electrons decreases, and thus the chemical potential (/j.f) is 
reduced, (c) When the chemical potential coincides with the 
mobility-edge, this model would predict that the MIT occurs, 
(d) The other end member of the homogeneity range of the 
spinel phase Lii+ ;c Ti2- : EO4(0 < x < |) is an empty band 
insulator (note that the chemical potential appears at the 
bottom of the allowed energy bands, and thus all electronic 
states are unoccupied), (e) The so-called quantum percola- 
tion threshold is reached when all electronic states have a 
localized character (regardless of the location of the chemical 
potential). 
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FIG. 2: Averaged numerical density of states for 
Lii+^T^-xCU systems with different doping concentrations, 
and the corresponding location of chemical potentials. The 
DOS of the undoped system corresponds to the thermody- 
namic limit, whereas plots for the doped systems are aver- 
aged density of states for 8 3 conventional unit cells using 10 
different realizations of disorder. 

threshold, one is guaranteed to find a maximally con- 
nected cluster, and (possibly) several isolated clusters. 
By definition, all electronic states associated with an iso- 
lated cluster are localized, while for the maximally con- 
nected cluster both extended and localized states will be 
found. Similar to Anderson localization 25 of disordered 
systems, we expect that the eigenstates of the maxi- 
mally connected cluster are localized in the band-edge 
and are separated from extended states at the middle of 
the band by the so-called mobility-edge^ 5 - By increas- 
ing disorder more and more states become localized and 
the mobility-edge moves toward the centre of the band. 
However, as the concentration of cations (producing the 
disorder) is increased the density of itinerant 3c? electrons 
is decreased. As long as the system's chemical potential 
(determined by both the density of 3d electrons and elec- 
tronic density of states of the disordered system) is above 
the mobility edge, the system remains metallic, whereas 
if it was below the mobility edge, the states in the im- 
mediate vicinity of the Fermi surface would be localized 
and the system would display insulating behaviour. The 
critical dopant concentration at which the chemical po- 
tential and mobility edge meet thus identifies the MIT. 
The diagrams in Fig. ^ summarize this process for the 
particular example of Lii +a; Ti2_ a ;04(for which x c ~ 0.15 
and xbi ~ 0.33). 

Example density of states curves, and the locations 
of the chemical potentials, are shown in Fig. [3 for differ- 
ent doping concentrations, averaged over systems with 8- 
cubed conventional unit cells (the largest studied). Note 
that the full of spectrum of eigenvalues, of both the max- 
imally connected cluster and all isolated clusters, have 



FIG. 3: Scaling of the relative localization length (A) of dif- 
ferent energy "bins" relative to that of the maximally con- 
nected cluster (Ao), for energies close to the mobility edge, 
for a doping concentration of x — 0.7. This ratio is plotted 
vs. the reciprocal of the number of conventional cells, and 
as the size of the lattices is increased from 2 3 conventional 
unit cells to 8 unit cells this data shows that the eigenstates 
with energies above -3.9 have delocalized behaviour. For this 
system we estimate that the (lower) mobility edge is located 
at -4.0 ±0.15. 

been evaluated since the latter contribute to the location 
of the chemical potential. 

The mobility edge has been estimated using the above- 
mentioned scaling method j2£ and Fig.|3depicts the appli- 
cation of this method to identify the location of mobility 
edge for the specific doping concentration of x = 0.7. 
To be concrete, one can estimate that for this dopant 
concentration the mobility edge, in units of the hopping 
energy t, is -4.0±0.15. 

We have combined all of our data in Fig. which 
shows the chemical potential and mobility edge that we 
would estimate for both the excess Li and doped Al sys- 
tems. For reference we have included the values of the 
band minimum for all dopings. Note that the Fermi 
level of the excess Li system crosses the mobility edge at 
roughly x c w 0.324, whereas for the doped Al system this 
crossing occurs at y c ~ 0.83. Clearly, these numbers are 
much larger than the experimental values {xmit ~ 0.15 
and dmit ~ 0.33). Thus, we have studied a system that 
should be very well represented, in a one-electron theory, 
by a quantum site percolation model, have determined 
the concentrations at which the Fermi levels cross the 
mobility edge, and find these values to be a factors of 2 
and greater than 3 larger than the experimental results. 
So, we believe that this shows that a one-electron model, 
that ignores electronic correlations, cannot explain the 
observed MITs. 
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FIG. 4: This plot displays the final numerical results of our 
study, and leads to the identification of predicted concentra- 
tions at which the MITs occur. These data are the the chem- 
ical potentials ((if) and mobility edge as a function of dop- 
ing for both Lii+ a; Ti2- : r04 and LiAl !/ Ti2-j / 04; for reference 
we have also plotted the minimum of the band of electronic 
states. The crossing of the chemical potential and mobility 
edge (denoted by open black circles) would indicate the po- 
sition of disorder-only induced metal-to-insulator transition, 
and these values are labelled x c and j/ c . 



Concluding, our results show that disorder-only mod- 
els of the MITs undergone by these systems substan- 
tially overestimate the critical concentrations of doped 
cations. To be specific, for Lii+^T^-zC^ our numerical 
result for little more than double the experimental 

value, and for LiAl a Ti2_ a 04 our numerical result is an 
even larger multiple. (We note that we have not elim- 
inated the possibility that polaronic effects give rise to 
this transition, although no experimental work points to 
their existence! 5 ) Thus, indirectly, this study supports 
the hypothesis that strong electronic correlations are im- 
portant for the MIT, and, possibly, also for the T c ~ 13K 
superconducting transition undergone by LiTi204. 

This past year— sample preparation advances have al- 
lowed for the growth of large stoichiometric LiTi2C>4 sin- 
gle crystals, and a remarkably sharp superconducting 
transition (8T C ~ 0.1 K) has been observed. We hope 
that these new samples stimulate further experimental 
research in this interesting class of materials. 
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